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1.  INTRODUCnON 

Determination  of  the  upper  atmospheric  wind  profile  between  80  and  200  km  by 
Observing  the  drift  of  vapour  trails  released  from  sounding  rockets  has  become 
well  established,  since  its  inception  in  1958  (Bedinger  et  al . , (ref . 1)) . Until 

recently  the  technique  was  restricted  to  twilight,  when  the  sky  was  dark  enough 
for  the  sunlit  trail  to  be  observed  from  the  ground.  In  order  to  extend  the 
period  of  observation,  and  to  provide  data  during  daylight  hours  when  many 
ionospheric  phenomena  are  more  evident,  an  instrument  to  discriminate  lithium 
against  the  daylight  sky  has  been  built.  A description  of  this  instrument  has 
been  given  by  Hind  and  Lloyd (ref . 2) , and  the  first  data  provided  by  it  are 
reported  in  reference  3. 

This  Note  presents  and  discusses  all  the  profiles  of  daytime  thermospheric 
winds  which  have  been  observed  since  the  initial  report. 

Section  2 discusses  the  data  analysis  methods,  and  Section  3 presents  the  data. 
In  Section  4 these  data  are  compared  with  profiles  determined  at  twilight,  and 
in  Section  5 the  calculated  redistribution  of  ionisation  expected  for  these  wind 
profiles  is  compared  with  observed  ground  based  observations  made  by  an  ionosonde 
and  magnetometer.  The  positions  of  these  instruments  with  respect  to  the  range- 
head,  along  with  those  of  the  daytime  scanner  observing  sites  and  typical  rocket 
trajectories  are  given  in  figure  1. 


2.  METHODS  OF  OBSERVATION  AND  DATA  ANALYSIS 

The  construction  and  operation  of  the  daylight  lithium  scanning  detector  has 
already  been  described  by  Hind  and  Lloyd (ref . 2) . The  instruments  operated 
satisfactorily  throughout  their  three  years  of  operation,  and  in  the  course  of 
their  operation  only  a few  minor  modifications  needed  to  be  introduced. 

The  main  modification  involved  automatically  recording  the  voltage  supply  to 
the  potentiometers  which  monitor  the  position  of  the  scanning  mirrors.  These 
reference  voltages  are  recorded  between  frames  during  reset,  and  are  processed 
in  the  same  manner  as  the  monitor  voltages  so  that  drifts  anywhere  in  the  record- 
ing system  are  taken  into  account.  It  was  found  that  the  reference  voltage  on 
the  final  processed  output  had  negligible  drift  during  the  trial  and  during  the 
angular  calibrations  (which  took  about  45  minutes),  but  that  differences  in  out- 
put equivalent  to  an  angular  change  of  up  to  two  degrees  could  occur  between  the 
trial  and  the  calibrations. 

Originally,  it  was  intended  to  use  the  stars  for  the  angular  calibrations, 
but  difficulties  in  distinguishing  the  star  data  from  spike  noise  (and  problems 
of  clouding  in  of  the  site),  necessitated  abandoning  this  method  in  favour  of 
using  the  sun.  The  azimuth  and  elevation  of  the  mount  on  successive  calibration 
points  were  changed,  so  that  a series  of  calibrations  were  made  with  the  image  of 
the  sun  lying  in  a different  position  in  the  field  of  view.  Thirteen  calibra- 
tion points  were  used,  chosen  to  lie  in  a matrix  evenly  covering  the  field  of 
view.  The  accuracy  of  the  calibration  procedure  was  checked  by  comparing  the 
actual  position  of  the  sun  with  the  sun's  position  calculated  from  the  calibra- 
tion parameters.  The  last  column  in  Appendix  I,  which  is  an  example  of  angular 
calibration  calculations,  shows  that  the  mean  difference  between  the  actual  and 
recomputed  sun's  |)ositions  is  about  0.3°.  Contributions  to  this  error  stem 
equally  from  the  signal  output,  the  scanner  mount  azimuth/elevation  readings, 
and  estimation  of  the  centre  of  the  sun  on  the  output  trace. 

The  accuracy  of  this  method  of  determining  the  wind  profile  has  been  checked 
against  the  proven  method  using  F24  cameras,  by  observing  lithium  trails  released 
by  two  Corclla  rockets  at  twilight.  A comparison  of  the  two  sets  of  data  also 
enabled  improvements  in  the  analysis  of  data  from  the  daytime  scanner  to  be  made. 
As  described  in  reference  2,  the  final  wind  profile  was  originally  determined  by 
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feeding  the  data  into  a function  minimisation  routine  using  a powerful  method 
developed  by  Powell.  However,  it  was  found  that  the  routine  could  converge  on 
subsidiary  minima,  instead  of  the  true  minimum  represented  by  the  correct  wind 
profile.  Also,  to  find  the  true  minimum  required  much  hand  manipulation,  and 
was  demanding  in  computer  time.  In  its  place  a more  suitable  technique  was 
evolved  which  yielded  consistent  profiles  in  a relatively  straight  forward 
manner. 

The  principle  of  the  method  was  to  assign  altitudes  to  observations  of  the 
lithium  trail  at  similar  times  from  the  two  sites  so  that  they  were  consistent. 

As  an  example,  consider  the  data  taken  from  one  frame  recorded  at  each  site  (a 
"frame"  is  a complete  set  of  scans  through  the  field  of  view,  as  for  example, 
in  figure  8).  The  data  consist  of  the  azimuth  and  elevation  of  points  along 
the  trial,  and  the  time  of  observation.  If  we  knew  the  altitude  of  each  data 
point  at  one  of  the  sites,  the  North  and  West  position  of  this  point  in  the 
trail  could  be  calculated.  Since  the  trajectory  of  the  rocket  (which  gives  the 
initial  position  of  the  trail)  is  known  from  radar,  it  would  enable  the  drift 
velocity  of  the  trail  to  be  calculated.  With  this  knowledge,  it  would  be 
possible  to  calculate  where  in  space  the  trail  was,  at  the  time  data  was  being 
taken  at  the  other  site.  From  this,  the  elevation  and  azimuth  of  the  trail 
from  that  site  could  be  calculated,  and  the  results  should  agree  with  the  observed 
data.  In  fact,  we  do  not  know  the  altitude.  What  is  done  is  to  assign 
different  sets  of  altitudes  to  the  data  at  both  sites,  and  see  which  set  gives 
consistent  agreement  between  the  reprojected  azimuth  and  elevation  (calculated 
using  the  method  explained  above)  and  the  observations.  Generally,  constraints 
can  be  put  on  altitude,  so  that  the  first  estimate  was  generally  within  5 to  10 
km.  These  constraints  are  the  rocket  apogee  (known  from  the  radar  trajectory), 
the  altitude  of  well  defined  kinks  in  the  trail  (which  could  be  triangulated 
directly),  the  lowest  altitude  below  which  lithium  was  rapidly  oxidised  (about 
85  km)  and  comparison  of  the  earliest  observations  with  the  radar  trajectory  of 
the  rocket.  An  example  of  this  last  technique  is  given  in  figure  2.  This  was 
an  observation  of  a downleg  release  from  a Skylark  rocket,  and  the  distance  of 
the  observing  sites  from  the  release  did  not  give  a large  subtended  angle,  i.e. 
poor  triangulation.  The  figure  shows  that  the  altitude  where  the  loop  forms 
(c.f.  figure  3)  can  be  assigned  within  10  km  at  worst. 

An  example  of  the  new  method  of  analysis  is  shown  for  the  difficult  Skylark 
1271  release  in  figure  3.  It  is  seen  that  the  highest  altitude  has  been 
assigned  an  altitude  10  km  too  high,  the  altitudes  were  correct  in  the  middle 
range,  but  further  down  were  too  high  again.  The  graph  also  shows  that  at  low 
altitudes  the  data  from  the  two  sites  are  inconsistent  with  each  other,  i.e.  one 
or  other  or  both  sets  of  data  are  wrong.  The  trail  was  in  fact  very  faint  at 
these  altitudes  and  it  was  hard  to  pick  it  out  from  the  noise.  Comparisons 
such  as  figure  3 have  the  added  advantage  of  acting  as  a guide  to  show  which 
data  are  incorrect. 

Figure  4 shows  the  wind  velocities  calculated  for  trial  Lorikeet  2013  for 
which  data  were  obtained  from  only  one  site  (the  other  site  was  clouded  in) . 

In  this  case  there  is  no  check  on  the  assigned  altitudes  by  reprojection  to  the 
other  site  and  comparison.  However,  the  velocity  hodograph  is  fairly  insensitive 
to  assigned  altitudes,  the  problem  is  to  assign  altitudes  to  this  hodograph. 

Even  when  data  are  obtained  from  both  sites,  it  is  preferable  to  quote  the  accuracy 
of  the  derived  wind  profile  as  a mean  error  in  altitude  rather  than  in  velocity. 

A typical  value  for  this  error  is  ± 2 km,  although  the  error  in  the  Skylark  1271 
was  larger  than  this,  because  of  the  problems  mentioned  earlier. 
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3.  DATA  AND  RESULTS 

Table  1 lists  the  rocket  firings,  and  associated  information,  from  which 
successful  daytime  releases  of  lithium  vapour  were  made.  The  trail  was  observed 
from  both  sites  in  all  cases,  except  for  Lorikeet  2013  for  which  good  data  were 
obtained  from  one  site  only. 

Figure  6 shows  photographs,  taken  during  the  trial,  of  the  lithium  trail  on 
the  monitor  C.R.O.  screen.  Figures  5,  7 and  8 give  examples  of  computer 
reconstruction  of  data  recorded  on  magnetic  tape.  As  can  be  seen  in  figure  6, 
at  early  times  a coarse  and  fast  scan  is  used  until  the  operator  is  sure  he  has 
the  trail  in  his  field  of  view.  The  scan  is  then  made  finer,  so  that  accurate 
data  on  the  position  of  the  trail  may  be  obtained. 

On  several  trials  the  lithium  cannister  continued  burning  over  apogee;  for 
example,  up  and  downleg  trails  can  be  seen  clearly  in  figure  5.  As  mentioned 
earlier,  the  apogee  point  is  in  such  cases  clearly  identifiable. 

Figure  9 shows  the  evolution  of  the  observed  trail  from  one  of  the  sites  for 
trial  Cockatoo  3010.  The  border  of  the  graph  is  drawn  to  enclose  the  field  of 
view  of  the  instrument.  It  is  seen  that  over  the  period  of  observation, 
generally  lasting  15  minutes  (by  which  time  all  but  the  dense  apogean  portion  of 
trail  has  grown  too  faint),  excellent  data  on  the  evolution  of  the  trail  were 
obtained.  Even  when  overall  drift  of  the  trail  required  changing  the  azimuth 
and  elevation  of  the  scanner  mount  during  the  trial  in  order  to  keep  the  trail 
in  the  field  of  view,  equally  good  data  were  obtained. 

Figures  10  through  22  give  the  final  wind  profiles  for  all  firings,  in  hodo- 
graph  form. 


4.  COMPARISON  OF  RESULTS  WITH  THEORY 

All  the  daytime  wind  profiles  showed  counter-clockwise  rotation  with  increas- 
ing altitude,  which  is  the  common  feature  of  all  twilight  vapour  releases 
observed  at  Woomera , with  the  exception  of  Lorikeet  2015  whose  hodograph  was 
predominantly  polarised  in  the  North-South  direction.  The  classical  Hough 
function  for  the  diurnal  (1,  -2)  mode  shows  just  such  a behaviour  at  30°  latitude, 
i.e.  the  East-West  component  vanishes,  and  suggests  that  on  this  occasion  this 
mode  was  completely  dominant.  However,  it  cannot  be  as  simple  as  that,  because 
the  observed  wind  profile  changes  from  Northward  to  Southward  at  135  km,  which 
the  (1,  -2)  mode  does  not. 

Forbes  and  Garret (ref .4)  have  calculated  the  solar  diurnal  tidal  structure  in 
the  thermosphere.  By  comparing  the  latitudinal  variation  of  the  wind  at 
different  altitudes  with  the  classical  Hough  modes,  they  found  that  the  (1,1) 
mode  was  dominant  up  to  about  110  km  altitude,  but  that  the  (1,  -2)  mode  became 

dominant  above  that  altitude.  This  is  consistent  with  the  observed  behaviour  of 

L2015;  also,  as  will  be  discussed  later,  the  presence  of  the  (1,1)  mode  is  in 
keeping  with  the  observed  data  from  the  other  firings. 

For  all  of  the  other  profiles,  which  do  not  show  such  an  idiosyncratic 
behaviour,  no  clear  identification  with  any  particular  mode  may  be  made.  The 
commonly  observed  anti-clockwise  rotation  of  the  wind  vector  with  increasing 
height,  with  a vertical  wavelength  greater  than  30  km  is  attributed  to  a combina- 
tion of  the  solar  semidiurnal  (2,2),  (2,4)  and,  to  a lesser  extent  (2,3)  modes. 

As  Hong  and  Lindzen(ref .5) , who  developed  a three  dimensional  model  of  the  solar 
semidiurnal  tide,  have  shown,  minor  changes  in  any  one  of  these  modes  can  make 
gross  differences  in  the  final  profile.  The  mode  strengths  in  turn  depend 
critically  on  the  mean  winds,  especially  below  ISO  km  altitude.  It  is  partly 
for  this  reason  that  it  is  difficult  to  compare  our  data  with  their  predictions. 

It  is  also  because  Hong  and  Lindzen's  model  at  30°  latitude  gives  very  low 
velocities  at  100  km  which  increase  without  much  change  in  phase  to  about 
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70  ms*'  at  120  km,  at  which  altitude  the  characteristic  counter-clockwise  rotation 
of  the  slowly  changing  velocity  vector  starts  to  occur.  The  bulk  of  our  results 
lie  in  the  range  100  to  120  km  altitude  and  already  show  the  developed  spiral  at 
these  altitudes.  Almost  certainly  this  behaviour  in  our  data  is  due  to  con- 
tribution from  the  solar  diurnal  (1,1)  mode  which,  unlike  the  (1,  -2)  mode  at 
these  altitudes,  is  not  plane  polarised. 

The  only  firings  for  which  we  have  data  at  higher  altitudes  is  SL1207  (figure 
13)  and  SL1271  (figure  19),  which  were  fired  at  the  same  time  and  season,  although 
two  years  separated  them.  In  view  of  the  variability  of  the  upper  atmospheric 
tides  and  other  wind  components,  there  is  a remarkable  similarity  between  the 
profiles.  Both  have  a southwa''d  velocity  at  100  km,  rotating  counter-clockwise 
and  becoming  Westward  at  130  km  altitude.  Interestingly  enough,  the  model  of 
Hong  and  Lindzen  also  predicts  that  the  maximum  Westward  velocity  at  130  km  will 
be  at  about  1100  hr.  The  profile  for  SL1271  continues  to  follow  the  model  above 
these  altitudes,  but  that  for  SL1207  appears  to  be  perturbed  in  the  North-South 
direction.  The  nature  of  this  perturbation  suggests  it  is  due  to  contributions 
from  the  solar  diurnal  (1,  -2)  mode. 

The  rockets  C3010,  L2011  and  L2013  were  fired  at  approximately  two  hour 
intervals.  A comparison  of  the  observed  wind  profiles  shows  them  to  be 
remarkably  similar  in  shape,  but  displaced  successively  in  altitude  ( figures  14, 
16,  18).  Because  the  phase  of  the  tidal  modes  descend  in  altitude  with  time, 
it  would  be  expected  that,  if  one  mode  dominates  a given  wind  profile,  an 
orderly  descent  in  the  profile  with  time  would  be  observed.  However,  the 
structure  of  these  three  profiles,  especially  the  cusp  which  appears  at  112  km 
on  C3010,  suggests  that  no  single  mode  dominates,  but  yet  the  three  profiles 
maintain  very  similar  shape  as  the  waves  descend.  Since  the  modes  which 
combine  to  form  the  profile  have  different  vertical  phase  velocity,  it  leads 
one  to  expect  a change  in  profile  with  time.  We  have  no  explanation  for  this 
contradiction  between  expected  behaviour  and  the  observations,  other  than  to 
suggest  it  may  be  a non-linear  "lock-in"  of  the  modes. 


5.  WIND  INDUCED  IONOSPHERIC  EFFECTS 

During  the  latter  part  of  the  period  over  which  the  firings  were  made,  an 
ionosonde  and  magnetometer  were  operational  at  the  B5  site.  Both  instruments 
gave  data  for  the  weeks  starting  22  February,  24  May  and  28  June,  1976.  There 
were  no  magnetic  storms  recorded  for  at  least  24  hours  before  any  of  the  firings; 
the  sum  of  the  3-hourly  Kp  magnetic  index  for  the  previous  24  hours  being  9,  12 
and  10  respectively.  However,  the  ionograms  showed  T.I.D.  (travelling  iono- 
spheric disturbance)  activity  for  24  February  1976.  This  was  especially  strong 
from  0800  to  0900  C.S.T.  Figure  23  is  a tracing  off  the  ionograms  for  this 
period.  Note  that,  contrary  to  what  frequently  occurs,  these  T.I.D.s  did  not 
give  rise  to  sporadic  E layers. 

A more  general  mechanism  for  the  formation  of  sporadic  E layers,  is  that  of 
a vertical  shear  in  the  horizontal  wind  (Whitehead(ref .6)) . There  is  substantial 
evidence  that  temperate  zone  sporadic  E layers  are  produced  by  metal  ions, 
deposited  by  meteors,  which  have  been  concentrated  into  layers  by  the  dynamo 
action  of  the  horizontal  neutral  wind  across  the  earth's  magnetic  field  (see, 
for  example,  Narcisi (ref . 7) ) . The  layers  form  where  the  vertical  ion  motion 
decreases  through  zero  with  increasing  altitude.  We  have  calculated  the 
vertical  ion  motion,  v^,  using  the  expression 
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where 


P 


r r 
y.  z 


u 


X’ 


u 


y 


ratio  of  collision-frequency  to  gyro- frequency  for  ions 
cosine  and  sine  of  angle  of  magnetic  dip 

horizontal  component  of  neutral  wind,  in  magnetic 
Eastward  and  Northward  directions. 


We  have  had  to  neglect  the  contribution  of  the  ambient  electric  field  as  it  is 
not  known.  However,  this  contribution  is  not  greater  than  10%. 

In  figure  15,  17,  20,  22,  are  plotted  the  vertical  ion  velocities  of  the 
trials  for  which  there  were  ionograms  recorded;  C3010,  L2011,  SL1271  and  L2015 
respectively.  All  the  profiles  had  a zero  crossing  with  increasing  altitude  in 
the  sense  required  to  produce  ion  convergence,  the  respective  altitudes  being  121, 
95,  123  and  134  km.  Despite  this  favorable  condition  for  sporadic  E,  with  the 
exception  of  two  weak  layers  (and  some  very  strange  layers  produced  by  the 
release  which  will  be  discussed  in  another  report),  no  sporadic  E was  produced. 

The  exceptions  were  a weak  cusp  Es  layer  at  120  km,  at  around  the  time  at  which 
C3010  was  fired,  and  a broad,  broken,  and  variable  layer  lying  between  120  and 
130  km  from  an  hour  before  the  time  of  launching  of  SL1271,  until  just  before 
launching  time.  The  heights  of  these  layers  agree  with  the  predictions  of  the 
wind  sliear  theory,  but  their  weakness,  and  the  absence  of  any  sporadic  E on  the 
otJier  occasions,  indicate  that  favorable  wind  shears  are  not  a sufficient 
condition  for  Es,  but  that  the  presence  of  metal  ions  in  the  appropriate  region 
is  also  necessary. 


6 . SUMMARY 

This  collection  of  data  on  wind  profiles  in  daytime,  shows  similar  structure 
and  variations  as  have  been  previously  observed  at  twilight.  Particular  points 
that  come  out  of  the  observations  were: 

(a)  In  one  case  the  velocity  profile  was  polarised  in  the  North-South 
direction.  This  behaviour  is  predicted  for  the  solar  diurnal  tide 
(1,  -2)  mode,  which  suggests  that  this  mode  was  the  dominant 
contribution  to  the  observed  wind  profile. 

(b)  In  most  of  the  other  cases  the  major  contribution  to  the  observed 
wind  profiles  appeared  to  be  the  solar  semidiurnal  tide  above  120  km, 
and  the  solar  diurnal  (1,  1)  mode  below  that  altitude. 

(c)  The  two  Skylark  trials  were  fired  at  the  same  time  of  day  and  season, 
albeit  in  different  years.  They  gave  data  with  remarkable  similarity, 
which  agreed  well  with  the  predicted  solar  semidiurnal  tide. 

(d)  Although,  in  every  case  examined,  the  form  of  the  wind  profile  was 
favorable  for  the  formation  of  sporadic  E,  on  only  two  occasions  was 
Es  observed,  and  then  only  very  weakly.  The  winds  are  therefore  only 
a necessary,  but  not  a sufficient  condition  for  Es  formation;  it 
being  often  the  case  that  the  metal  ion  concentration  is  not  great 
enough  to  give  an  observed  layer. 
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TABLE  1.  DAYTIME  LITHIUM  RELEASES  AT  WOOMERA 


Vehicle 

Date 

C.S.T. 

U.T. 

Cockatoo  2004 

21  June  1973 

08  50  00 

23  20  00 

Cockatoo  2001 

19  September  1973 

14  30  00 

05  00  00 

Cockatoo  3004 

19  March  1974 

12  28  00 

02  58  00 

Skylark  1207 

23  April  1974 

10  33  00 

01  03  00 

Cockatoo  3010 

24  February  1976 

07  00  00 

21  30  00 

Lorikeet  2011 

24  February  1976 

09  00  00 

23  30  00 

Lorikeet  2013 

24  February  1976 

12  00  00 

02  30  00 

Skylark  1271 

27  May  1976 

10  00  00 

00  30  00 

Lorikeet  2015 

30  June  1976 

10  00  00 

00  30  00 

APPENDIX  I,  ANGULAR  CALIBRATION  OF  DAYTIME  LITHIUM  SCANNER  USING  THE  SUN 
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Comparison  of  observed  data  with  reprojected  data  from  other 
site  using  method  of  assigned  altitudes 


Figure  4.  MTind  hodograph  calculated  from  data  from  site  B8 
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Note:  Release  of  lithium  continued  over  apogee.  These 

computer  plots  were  generated  from  magnetic  tape  data. 


Figure  5.  Observations  of  lithium  trail  from  Cockatoo  3004 
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These  photographs  were  taken  off  the  monitor  CRO.  The  perturbation  at  the  bottom  of  frame  14  is 
due  to  cloud. 
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Figure  7.  Observations  of  litnium  trail  from  Cockatoo  3010 
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Figure  8 


Observation  of  lithium  trail  from  Lorikeet  2011 
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For  determining  the  Chebychev  coefficients  the  profile  and  the 
data  were  separated  into  two  height  bands.  The  connecting 
points  are  plotted  as  circles. 


Figure  13.  Wind  hodograph  for  Skylark  1207 
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Note:  Altitudes  could  not  be  assigned  as  data  were 
obtained  only  at  one  site. 


Figure  18.  Wind  hodograph  for  Lorikeet  2013 
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Figure  22.  Vertical  ion  motion  calculated  for  Lorikeet  2015 
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